Phase fluctuations versus Gaussian fluctuations in optimaUy-doped YBa2Cu307. 
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We analyze recent torque measurements of the magnetization A'ld vs. magnetic field H in optimally 
doped YBa2Cu307-y (OPT YBCO) to argue against a recent proposal by Rey et al. that the 
magnetization results above Tc are consistent with Gaussian fluctuations. We find that, despite its 
strong interlayer coupling, OPT YBCO displays an anomalous non-monotonic dependence of Md 
on H which represents direct evidence for the locking of the pair wavefunction phase On at Tc and 
the subsequent unlocking by a relatively weak H . These unusual features characterize the unusual 
nature of the transition to the Meissner state in cuprate superconductors. They are absent in low-Tc 
superconductors to our knowledge. We also stress the importance of the vortex liquid state, as well 
as the profiles of the melting field Hm{T) and the upper critical field curve Hc2{T) in the T-H plane. 
Contrary to the claims of Rey et al., we show that the curves of the magnetization and the Nernst 
signal illustrate the inaccessibility of the Hc2 line near Tc- The prediction of the Hc2 line by Rey et 
al. is shown to be invalid in OPT YBCO. 



I. INTRODUCTION 

In the mean-field Gaussian treatment of fluctuations 
(valid in low- Tc superconductors), the pair-wavefunction 
amplitude \'^\ vanishes at the critical temperature Tc. 
Above Tc, amplitude fluctuations about the equilibrium 
point I^Pl = may be regarded (in Schmid's elegant de- 
piction [l|) as droplets of condensate of radius ^gl, the 
Ginzburg-Landau (GL) coherence length. In the compet- 
ing phase-disordering scenario, j^*! remains flnite above 
Tc. The collapse of the Meissner effect above Tc is caused 
instead by the vanishing of phase rigidity. Above Tc, fluc- 
tuations are primarily of the phase 6* of proceeding by 
phase slips caused by the motion of spontaneous vortices. 
In zero magnetic fleld, the net vorticity in the sample is 
zero, so the populations of "up" and "down" vortices are 
equal. In underdoped cuprates, there is strong evidence 
from Nernst 0, magnetization and other experi- 

ments in support of the phase-disordering scenario. 

Initially, it seemed to us that optimally-doped 
YBazCugOr-j, (OPT YBCO), which has the largest 
interlayer coupling energy (and lowest electronic 
anisotropy) among the known cuprates, would be most 
amenable to standard mean-field (MF) treatment, i.e. 
its fluctuations above Tc are strictly Gaussian. In- 
deed, in the early 90s, several groups applied the con- 
ventional Maki- Thompson Aslamazov-Larkin theory to 
analyze fluctuation conductivity above Tc 6]. Sub- 
sequently, these mean-field "fits" , seemingly reason- 
able in OPT YBCO, were found to be woefully inade- 
quate in underdoped YBCO and in most other cuprates. 
The steady accumulation of evidence from Nernst 
and torque experiments favoring the phase-disordering 
mechanism in underdoped hole-doped cuprates and in 
OPT Bi2Sr2CaCu208+j; (Bi 2212) has prompted a re- 
assessment of the case for OPT YBCO. The recent torque 
magnetometry results in Ref. Q have persuaded us that 
OPT YBCO is, in fact, much closer to the other hole- 



doped cuprates. The collapse of its Meissner state at Tc 
is also caused by vanishing of phase rigidity. 

In a Comment, Rey et al. have fitted the magneti- 
zation curves for OPT YBCO in Ref. [Sj to their model 
and claimed that the diamagnetic signal is consistent 
with MF Gaussian fluctuations in a layered superconduc- 
tor. Furthermore, they have extended their calculation to 
temperataures T < Tc to infer an upper critical field Hc2 
that rises linearly with reduced temperature (1 — T /Tc) 
with a slope of -3 T/K. 

Here, we show that the fluctuation signal above Tc is 
just the tip of the iceberg. Because of strong interlayer 
coupling in OPT YBCO, it is necessary to go below Tc to 
uncover the close similarities with other cuprates. In Sec. 
|TT]we describe the non-monotonic variation of the magne- 
tization curve just below Tc, previously noted in torque 
measurements on OPT Bi 2212 This characteristic 
feature - inherent in the phase-disordering scenerio - is 
common to the cuprates investigated to date. In Sec. 
mil we discuss the upper critical fleld Hc2 for T close Tc 
in OPT YBCO. The striking inability of high-resolution 
experiments to detect the Hc2 fine, a corner stone of the 
mean-field Gaussian picture, is also a characteristic fea- 
ture of the phase disordering mechanism. Finally, in Sec. 
IIVI we discuss the fits to the magnetization above Tc- We 
show that the conclusions of Rey et al. Q are not valid. 



II. NON-MONOTONIC MAGNETIZATION 

The curves of the diamagnetic component of the mag- 
netization Md in OPT YBCO (Ref. [5|]) are reproduced 
in Fig [TJ At temperature T above ^100 K, Md is ini- 
tially linear in up to 6 T. As we approach Tc, the 
linear-response regime becomes confined to progressively 
smaller field ranges (for e.g. \H\ < 2 T at 95 K). 

At Tc = 92.5 K (the present analysis shows that Tc is 
slightly higher than the nominal value in Ref. [5]), the 
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FIG. 1: (color online) The field dependence of the diamag- 
netic component of the magnetization Md in optimally doped 
YBa2Cu307_y at selected temperatures T. At Tc — 92.5, a 
weak step-increase in \Md \ signals the onset of full flux expul- 
sion. Below Tc, a field H > Hcusp causes a sharp decrease in 
the screening current, but at larger H, \Md\ resumes its in- 
crease at a rate similar to that above Tc. This non-monotonic 
pattern signals the field-induced decoupling of 9n between ad- 
jacent bilayers below Tc- We argue that this is the defin- 
ing magnetization feature that characterizes the transition in 
cuprates. no is the vacuum permeability. The arrows (la- 
belled as //c2?) are Hc2 values at 89, 90 and 91 K predicted 



inRef. 0]. (Adapted from Ref. d.) 



only significant change is a step-increase in the magni- 
tude \Md\ near H ^ 0+ . Above 0.5 T, Md{H) is strik- 
ingly similar to the curves above Tc apart from vertical 
scale. We focus on the region in Fig. [T] bounded by the 
curves at Tc and 90 K (the curve at 91 K is represen- 
tative). At H — 0^, the step-increase in \Md\ signals 
the onset of full flux expulsion. However, a weak field 
interrupts this steep rise to produce a sharp cusp at the 
field Hcusp {Hcusp is slightly higher than the lower critical 
field Hci because of surface pinning of vortices). When 
H exceeds Hcusp, \Md\ falls rapidly to a broad minimum 
near 1 T, but subsequently rises to even larger values. 
This non-monotonic profile in Md, appearing just below 
Tc, is a defining hallmark of the transition in cuprates. In 
this interval (90-92.5 K), Md vs. H is reversible except 
in the vicinity of Hcusp, where surface-barrier pinning of 
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FIG. 2: (color online) The magnetization curves below Tc in 
NbSe2. At each T, Md{H) decrease monotonically to reach 
zero at Hc2(T). Following vortex entry at Hd (not resolved 
here), the diamagnetic screening currents are steadily weak- 
ened as H suppresses the pair amplitude. The inset plots the 
inferred Hc2 vs. T. (Adapted from Ref. Q|) 



vortices causes a slight hysteresis. 

This unsual non-monotonic profile - absent in low-Tc 
superconductors to our knowledge - implies the follow- 
ing picture. The collapse of the Meissner state at Hcusp 
leads to a steep decrease in \Md\ as vortices enter the 
sample. This is a consequence of field suppression of 
the nascent inter-bilayer phase coupling. (The two Cu02 
layers bracketing the Y ions constitute a bilayer. We are 
concerned with only the coupling between adjacent bilay- 
ers which are separated by a spacing d— 11.8 A. We men- 
tion the intra-bilayer coupling below.) In low-Tc super- 
conductors, the decrease invariably continues monotoni- 
cally to zero as H Hc2{T) (the upper critical field at 
temperature T). For reference, we show in Fig. [5] curves 
of Md vs. H in the layered superconductor NbSe2 d]. 
By contrast, the non-monotonicity in \Md\ implies that 
the diamagnetic current in OPT YBCO turns around 
and grows stronger with _ff at a rate closely similar to 
that above Tc- Remarkably, if we hide the field region 
around Hcusp, the profiles of Md vs. H below Tc resem- 
ble those above Tc, apart from a different vertical scale. 
This implies that the transition at Tc at which flux ex- 
pulsion appears is a weak-field phenomenon. Beyond this 
weak-field regime, there is no hint that a transition has 
occured. 

Returning to Fig. [U we interpret the unusual non- 
monotonic pattern (now seen in Bi 2201, Bi 2212, LSCO 
and YBCO) [1-[H as reflecting the rapid growth of c- 
axis phase rigidity below Tc (in zero or weak H) and 
its subsequent destruction at Hcusp- In a finite interval 
above Tc, the pair wavefunction |^'n|e'^" in bilayer n has 
a finite average amplitude (|5'„|). The in-plane phase- 
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phase correlation length Cq in each bilayer, given by 

is long enough that local diamagnetic currents can be 
detected by torque magnetometry. However, the c-axis 
correlation length ^ d, so On is uncorrelated between 
adjacent bilayers. Hence Md reflects the diamagnetic re- 
sponse of 2D supercurrents which are observable to very 
intense H . In Fig. [Tl this 2D diamagnetic response is 
represented in the field profile of Md at 94 K. 

Below Tc, diverges (in zero H) to lock 9n across all 
bilayers. In weak H, the 3D phase rigidity produces full 
expulsion. However, in the 2-K interval 90-92.5 K, a weak 
field H ~ Hcusp suffices to destroy the c-axis phase stiff- 
ness. The system then reverts back to the diamagnetic 
response of 2D uncorrelated condensates. Consequently, 
as H increases further, the 2D diamagnetic response con- 
tinues to increase, mimicking the profile at 94 K. At 14 T, 
the 2D response leads to a value for \Md\ that well exceeds 
its value at Hcusp (see 91-K curve). The weak locking of 
9n across bilayers, and the field-induced unlocking, ac- 
count for the non-monotonicity of Md as well as the sim- 
ilarities of the high-field portions across Tc in a physically 
reasonable way. The juxtaposition of an extremely large 
pairing energy scale (d-wave gap amplitude A ■~40 mV) 
and a weak c-axis phase stiffness leads to this very un- 
usual phase locking- unlocking scenario which seems per- 
vasive in the hole-doped cuprates. We argue that this 
situation cannot be treated by applying the mean-field 
GL approach [l| to the Lawrence-Doniach (LD) Hamil- 
tonian 

III. CLOSURE OF THE Hc2 CURVE 

Another strong argument against the mean-field de- 
scription is obtained from the qualitative features of the 
phase diagram in the T-H plane (Fig. [3|) . A cornerstone 
of the mean-field Gaussian description of a type-II super- 
conductor is the well-known profile of the upper critical 
field curve -ffc2(T) (Fig. Hi). In the T-H plane, the Hc2 
curve is a sharply defined boundary that separates the re- 
gion with finite pair amplitude 1 4* | from the normal state 
with \^\ =0. The Gaussian fluctuations are fluctuations 
of the amplitude around the |^| = state above Hc2{T). 
In low-Tc superconductors (for T 0.5 T^), Hc2{T) de- 
creases linearly in the reduced temperature t = 1—T/Tc 
to terminate at T^. This proflle is shown for NbSe2 in 
the inset of Fig. [2] To stress this point, we refer to the 
termination as "closure" of the curve of Hc2- 

The closure of Hc2 reflects an anomalous aspect of the 
BCS transition at Tc. As T is elevated above Tc in zero 
i?, the pair condensate vanishes before it loses its phase 
rigidity. This implies that both the superfluidity (which 
depends on flnite phase rigidity) and the physical entities 
that manifest long-range phase coherence (the Cooper 
pairs) vanish at the same temperature T^. (This is akin 
to having the local moments in an antiferromagnet vanish 
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FIG. 3: (color online) Profiles of the melting field HmiT) and 
upper critical field Hc2{T) in low-Tc and cuprate superconduc- 
tors. Panel (a) shows the phase diagram for a conventional 
type-II superconductor in the T-H plane. The vortex liquid 
phase is wedged between the curves of Hm{T) and Hc2{T), 
both of which terminate at Tc as H ^ 0. In the hole-doped 
cuprates [Panel (b)], the vortex liquid state dominates the 
phase diagram. As iif — >■ 0, HmiT) terminates at the observed 
Tc whereas Hc2{T) terminates at a higher temperature Tmf- 
Panel (c) : The contour plots of the magnitude of the Nernst 
signal Cy in OPT YBCO in the T-H plane. Values of the 10 
contours are displayed on the left column. The melting field 
Hm{T) separates the vortex solid (black region) from the vor- 
tex liquid state. The dashed line is the Hc2 line predicted in 
Ref. [3. (Adapted from Ref. [H.) 



at the Neel temperature.) In an applied field H, phase 
rigidity is maintained if the vortices are in the solid state 
(and pinned). As H approaches Hc2, the phase rigidity 
abruptly vanishes at the melting field Hm which marks 
the transition to the vortex-liquid state (the analog of the 
paramagnetic state in the antiferromagnet). The vortex- 
liquid state is wedged between the curves of Hm{T) and 
Hc2{T) (Fig. 13^). As sketched, the anomalous feature 
in the BCS scenario is that, as _ff ^> 0, the vortex-liquid 
region vanishes. 

The competing phase-disordering scenario has a qualti- 
tatively different phase diagram (Fig. IJb)- The vortex 
liquid state now occupies a much larger fraction of the 
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FIG. 4: (color online) Plots of the Nernst signal ey = Ey/\\JT\ 
vs. H in OPT YBCO with Tc = 92 K at selected T {Ey is the 
Nernst i5- field and — VT the applied temperature gradient). 
At the melting field Hm, Cy rises nearly vertically reflecting 
the sharp increase in vortex velocity. The contour map (Fig. 
O;) was derived from these curves (and additional curves not 
plotted). The "up" arrows denote the values for Hc2 at 88 
and 90 K predicted in Ref. 0]. (Adapated from Ref. 1(|.) 



region in which ^ 0. Significantly, the melting curve 
HraiT) intercepts the H = axis at a temperature Tm{0) 
lower than the termination point Tmf of the Hc2 curve, 
as sketched in Panel (b). If we increase T along the 
H — axis, phase rigidity is lost at Tm(0) long before 
we reach Tmf- Since both the Meissner effect and the 
zero-resistivity state are crucially dependent on having 
long-range phase stiffness, the Tc commonly observed by 
these techniques is identified with Tm(0). 

Above Tc, we have a vortex liquid whose existence may 
be detected using the Nernst effect and torque magne- 
tometry (resistivity is ineffectual). The onset tempera- 
tures Tonset for the vortex Nernst and diamagnetic sig- 
nals 0, Q are lower bounds for Tmf- Clearly, if the 
interval \Tc,Tmf] is large (underdoped cuprates), the 
curve of Hc2{T) is nearly T independent in the inter- 
val > T > Tc, it attains closure only at Tmf- (The 
T-independent part has been established in underdoped, 
single-layer Bi 2201 where Hc2 (~50 T) is nearly acces- 
sible with a 45-T magnet 0.) Hence, if we focus on 
experiments close to Tc, the only curve that intersects 
the H = axis is the melting curve. Even at Tc, Hc2 is 
inaccessible except for Bi 2201. The absence of closure 



for Hc2{T) implies that the vortex-liquid state extends 
above Tc to Tmf, and survives to H Hc2- 

A striking empirical fact in the hole-doped cuprates is 
that experiments (over a 25-year period) have never ob- 
served an Hc2 curve that decreases linearly with t to ter- 
minate at Tc (not counting early flux-flow resistivity ex- 
periments that mis-identified Hm for Hc2)- The absence 
of the Hc2 curve is incompatible with the Gaussian pic- 
ture, but anticipated in the phase-disordering scenario. 

We now turn to results in OPT YBCO. Figure [3J: dis- 
plays the contour plots of the Nernst signal ey in another 
OPT YBCO crystal with identical Tc an d closely com- 
parable quality (adapted from Ref. [li|). The Nernst 
signal is defined as ey = Ey/\yT\, where Ey is the trans- 
verse electric field produced by the velocity of vortices 
diffusing in the applied thermal gradient — VT (see Ref. 
0). In the vortex solid {H < HmiT)), ey is rigorously 
zero. Just above Hm{T), ey rises very steeply reflecting 
the steep increase in vortex velocity in the liquid state. 
In comparison with similar contour plots for underdoped 
cuprates (as well as OPT Bi 2212), the interval [Tc, Tmf] 
in OPT YBCO is smaller. Also, the curve Hm{T) rises 
from Tc with a much steeper slope. Nevertheless, the 
contour features are qualitatively similar. The contours 
are very nearly vertical at Tc, implying that ey hardly 
changes with H . Even exactly at Tc, there is no field 
scale above which the system can be described as being 
in the "normal state" with \^\ = 0. In cuprates, we 
find it helpful to regard the segment of the H — axis 
between Tc and Tmf as the continuation of the H„i{T) 
curve. The vortex liquid, confined between the curves 
of Hm and Hc2 below Tc, expands to occupy the entire 
region below Hc2 above Tc- It plays a dominant role in 
thermodynamic measurements. 

In Ref. [Tl , the i?c2 line is predicted to increase linearly 
with t from the point {T — Tc, H — 0) with a slope of 
-3 T/K. We have drawn their Hc2 line as a dashed line 
(labelled as Hc21) in Fig. [3J:. Clearly, the predicted Hc2 
line cuts across the contours in an arbitrary way that 
bears no resemblance to experiment. 

To explain this better, we show in Fi g. HI the profiles 
of ey vs. H at selected fixed T near Tc [lOI- Just below 
Tc (curves at 90 and 88 K), ey rises nearly vertically 
when H exceeds Hm, reflecting the sharp increase in the 
vortex velocity in the vortex liquid in response to the 
gradient — VT . The magnitude of ey, which remains 
large up to 14 T, implies that remains finite. There 
is no experimental feature (change in slope, e.g.) that 
signals I^E"! — >■ at the values of Hc2 predicted by Rey et 
al. (shown as "up" arrows). 

The same difficulty exists for the plots of Md vs. H in 
Fig. [Tl where the predicted Hc2 values are indicated by 
the arrows. Contrary to the claim in Ref. .l}, the curves 
of Md vary smoothly through the predicted field values 
at 89, 90 and 91 K. Again, there is no feature reflect- 
ing either a transition or crossover. Given that \Md\ is a 
measure of the supercurrent density Js surrounding each 
vortex in the sample (even in the vortex liquid state) , we 
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FIG. 5: (color online) Plots of the susceptibility x = M/H 
vs. ^.qH in OPT YBCO at selected T (T^ = 92.5 K). The 
dashed horizontal lines indicate the linear-response region for 
curves near Tc- 



would expect \Md\ to be considerably larger when 1 5*1 is 
finite [H < Hc2) compared with the Gaussian fluctua- 
tion regime above Hc2- Instead, \Md\ hardly varies as H 
crosses the predicted values. One confronts the very awk- 
ward problem of explaining why | Md \ retains nearly the 
same value below and above "iJc2"- Even worse, close 
to Tc (curve 91 K), \Md\ actually increases when H ex- 
ceeds 3 T (this reflects the decoupling of the bilayers as 
discussed in Sec. HI)) . 



IV. FLUCTUATIONS ABOVE 

We comment on the fits above reported in Ref. Q. 
Figure [5] shows plots of the susceptibility x — M/H. 
Above 100 K, x is -ff-independent over a large field 
region. However, as T falls below 100 K, the linear- 
response is confined to progressively smaller field ranges, 
as mentioned. Nonetheless, as T — > T+, the linear- 
response value of X (dashed lines in Fig. [5]) does not 
diverge. Instead, the appearance of the Meissner effect 
(and its subsequent suppression by H) occur at much 
lower field scales as discussed. Below the hallmark 
non-monotonicity of Md vs. H described above is harder 
to see in plots of x vs. H because of the strong field 
variation caused by dividing a relatively fiat profile by 



FIG. 6: (color online) The T dependence of the susceptibility 
in OPT YBCO for selected values of H. Above Tc (92.5 K), x 
is H independent below 2 T (linear response region). Below 
Tc, X becomes strongly H dependent at low H because of the 
field suppression of phase locking in the vicinity of Hcusp- At 
large H (curves at 6 and 14 T), x varies smoothly through 
Tc. The dashed line is the mean-field (MF) fit to Eq. [5] with 
■n = 3.78 X 10"^ and Bld = 0.11. 



H. 

Moving to the fiuctuation regime above Tc, we show 
in Fig. [S] the plots of x vs. T at selected values of H. 
Above Tc, all values of x below 2 T collapse to a single 
curve, which defines the linear-response fiuctuation sus- 
ceptibility that can be directly compared with mean-field 
theories P, H, Q . For illustration, we show as a dashed 



curve the best fit to the linear-response MF expression 



xiT) 



eB 



(2) 



LD 



with r/ an adjustable numerical factor and e = [T —Tc) /Tc 
(the dashed curve has T] = 3.78 x 10""^). For the LD 
parameter, we used the value proposed by Rey et al. Q 
[Bld = 0.11). 

While the fit can account for the overall magnitude of 
X above 100 K with reasonable parameters, as shown by 
Rey et al., the functional form in Eq. [5] does not describe 
the data trend all that well. Between 94 and 98 K, it 
underestimates |x| by as much as 15%, and overestimates 
its value at 92.5 K. A similar pattern of overshooting and 
undershooting is also evident in Fig. 1(a) of Rey et al.. 
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At large H (curves at 6 and 14 T), x varies smoothly 
through Tc without detectable change in slope. This re- 
flects the continuity of the diamagnetic response in strong 
fields described above. Significantly, at low H, the linear- 
response segment of x suddenly becomes strongly H- 
dependent even in very weak H. As we discussed, this re- 
flects the non-monotonicity that appears above the field 
Hcusp- These features cannot be described by the Gaus- 
sian approach of Rey et al. Q . 



V. SUMMARY 

In attempting to understand the fluctuation signals of 
the magnetization in OPT YBCO, it is essential to view 
the results above and below Tc- We show that the non- 
monotonic curves of Aid vs. H just below Tc are con- 
sistent with the loss of long-range phase coherence in a 
layered superconductor with extremely large pair-binding 
energy within each layer, but a c-axis coupling (between 
bilayers) that is suppressed by a few Teslas close to Tc- 
These features are incompatible with the Gaussian pic- 
ture. In addition, when discussing the thermodynamics 
of hole-doped cuprates, we argue that it is vital to rec- 
ognize the elephant in the room, namely the vortex liq- 
uid above and below Tc- The dominant presence of the 
vortex liquid alters qualitatively the profiles of the Hc2 
curve, pushing its closure to a temperature higher than 



Tc- As a result, the corner-stone feature of the Gaussian 
approach, namely an Hc2 line that terminates at Tc, is 
absent in OPT YBCO (and other hole-doped cuprates). 
This contradicts a prediction of Ref. Q . 

In the phase-disordering picture, the d-wave pairing 
gap A must persist high above Tc- Spectroscopic exper- 
iments are increasingly able to distinguish between the 
gap in the vicinity of the node from the much larger 
antinodal gap. The persistence of the nodal A above 
Tc has now been reported in scanning tunneling mi- 
croscopy |12| and photoemission experiments 13, 14] on 
OPT Bi 2212, and in c-axis infrared reflectivity exper- 
iments on YBCO [l^. The vanishing of the Meissner 
effect above Tc reflects the collapse of the inter-bilayer 
phase-coupling. Since the intra-bilayer phase coupling is 
much stronger, one might expect that the concomitant 
phase rigidity can be observed experimentally above Tc- 
Recently, this was detected by Dubroka et al. [l^ as a 
Josephson plasma resonance that persists above Tc over 
a broad doping range in YBCO (as high as 180 K in the 
underdoped regime). 

We have benefitted from valuable discussions with 
Zlatko Tesanovic, Oskar Vafek, Steve Kivelson, Ashvin 
Vishwanath, Sri Raghu, Ali Yazdani, Christian Bernhard 
and P. W. Anderson. Support from the U.S. National 
Science Foundation under MRSEC grant DMR 0819860 
is gratefully acknowledged. 
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